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We constructed the one-particle spectral functions (diagonal and off-diagonal) which reproduce 
BCS for weak coupling and which take into account the effect of correlations on superconductivity 
in the attractive Hubbard model. The diagonal spectral function is composed of three peaks and 
the off-diagonal one is composed of two peaks. This ansatz satisfies the sum rules for the first 
six moments. Our solutions are valid for intermediate coupling, i.e., for U/t » —4.0. Our set of 
analytical equations for the unknown variables is self-consistent and has been solved numerically in 
lowest order of the order parameter. As a result, we obtain that the presence of the third band, 
or upper Hubbard band, strongly renormalizes the two lower bands, making that the energy gap be 
k-dependent while the order parameter is pure s-wave. This shows that the order parameter and 
, the gap are two different quantities. 
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The study of correlations has been renewed again by the discovery of high T c superconducting oxide materials 
(HTSC), jjj since these materials exhibit a short coherence length, £, and a very large penetration depth, A, thus 
behaving like extreme type-II superconductors, i.e. k = A/£ ^> 1. 
-H , We will use the on-site attractive Hubbard Hamiltonian put forward by Micnas et al |2j as a phenomenological 
model for describing the HTSC. Previous authors have used this model to study the bismuthate superconductors. ^ 
The importance of the Hubbard model has been recognized by Denteneer et al Q who say that if the phase diagram 
of the Hubbard model were fully understood, it might form the basis of understanding correlated electrons as much 
as the Ising model did for understanding critical phenomena. Schneider et al ^| have applied this model to explain 
universal properties of several families of these compounds, like the relation between transition temperature, magnetic 
\Q ' penetration depth and gap, at zero temperature. It is likely that understanding this will provide insight into the effect 
ON ■ of correlations on measurable quantities as shown by Singer et al || . 

We will use the exact relations of Nolting for the one-particle diagonal spectral function together with the exact 
relations for the one-particle off-diagonal spectral function, i.e., the anomalous Green's function, to study systems 
with broken symmetry. || We will start from the dynamical equations constructing one-particle Green's functions 
1 , I which reduce to the BCS solution in the case of weak coupling. Our solutions, even in the case of BCS, give in a 
1 ■ natural way the Hartree shift to the chemical potential. 
q The Hamiltonian we study is the following 
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where c\ a (ci G ) are creation (annihilation) electron operators with spin a, = c\ a Cic t is a hopping matrix element 
between the nearest sites i and j, and U is the onsite interaction, p, is the chemical potential and we work in the 
grand canonical ensemble. In the present study we consider an attractive interaction, U < 0. The Hamiltonian of 
Eq. ([!]) has been studied in detail by Micnas et al 

Nolting B and Kalashnikov and Fradkin || have derived exact relations for the one-particle spectral functions, 
A(k, ui) and B(k,uj), which are given by: 



M (0) (k) = / A{k,uj)duj = 1 , (2) 
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B(k,a;)dw = 



w5(k,w)da; = A Q (T) 
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where e(k) = — 2t(cos(k x a) + cos{k y a)), with a the lattice constant. A (T), the off-diagonal long range order, and p, 
the carrier concentration per spin, n/2, have to be calculated self-consistently. A(\t,u>) and S(k, ui) are the diagonal 
and off-diagonal one-particle spectral functions, respectively. They are defined by: 



A(k,w) 



— lim 7m[G(k, w + i5)] 
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The parameter B & , the Nolting's approximation, 
from 



BCk,u) = lim Im[F(k,u + i6)} . (8) 

is given in a self-consistent way. /? and A Q (T) are calculated 
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The main ingredient in any correct description of a physical system is the availability of the one-particle Green's 
function. Of course, in a many-body problem we cannot have a closed expression. In order to show the way, we first 
discuss the 2D BCS case in the framework of the moment approach and then we study the case of superconductivity 
beyond BCS, i.e., by including higher order moments. Here one must note that due to the Mermin- Wagner theorem 
no phase-transition is expected in 2D in systems with a continuous symmetry, and the formalism here is too simple 
to describe a Kosterlitz-Thouless phase-transition. Nonetheless, it is observed that the formalism does give a phase- 
transition. 



A. BCS case 



Let us postulate as our one-particle spectral functions the following ones: 

i4(k,w) = ai(k)S(oj - E k ) + a 2 (k)6(oj + E k ) , (10) 
B{k, u) = v/ai(k)a 2 (k) [5(u - E k ) - 6(u + E k )\ , (11) 

By using this Ansatz in the moment equations, we obtain that: 



£k = V(tk-V + pUr + Al(T) ; ttl (k) = -^l+ k ^T" j ; a 2 (k) = l-ai(k) . (12) 

Eqs.(^2|) are clearly the BCS solutions with the Hartree shift, pU , included in a natural way. The quasi-particle 
excitation spectrum, E kl has been obtained from the particular form of the off-diagonal single-particle spectral 
function, i.e., the square root dependence. This square root dependence will be kept when we include correlations, 
i.e., when we include higher order moments, i.e., a 2 and 03. This we will do next. 



B. Correlations on the BCS solution. 



Now, we will go beyond the BCS solution by including correlations in higher order. In this case, the physical picture 
is that superconductivity is similar to BCS in the sense that we will have quasi-particle spectra with a gap around 
the chemical potential. 

In order to keep the picture simple, we will assume that the chemical potential, p,, is close to the bottom of the 
free band. Then, the superconducting gap opens up around the chemical potential and there is an upper Hubbard 
band which remains almost similar to the Hubbard band when the order parameter is zero, i.e., when there is no 
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symmetry breaking. In this case, we postulate that the diagonal one-particle spectral density is composed of three 
poles as follows: 



A(k, to) = ai(k)S(u - Oi(k)) + a 2 (k)S(u - fi 2 (k)) + a 3 (k)S(uj - 3 (k)) 



(13) 



In Eq. (|lj), the first two poles represent the behavior around the chemical potential and the third pole is due to the 
influence of the upper Hubbard band, which in the case of the atomic limit, t — 0, is the upper Hubbard band or the 
band of single occupied states. The lower band, now split in two, is the band of doubly occupied states for the case 
of almost atomic limit ||. The off-diagonal (anomalous) spectral function has the same form given as in the BCS 
case. Of course, our parameters cty(k) and Oj(k) (j = 1,2,3), have to be calculated self-consistently (f2i(k) = +.E k ; 
fi 2 (k) = — -Ek! ^3(k) = fl k )- Now, the parameter B s which appears in a3(k) is generally k-dependent. In the 
spherical approximation , it is k-independent and is calculated from the following equation: 
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where (3 = l/(fcsT), T the absolute temperature, ks the Boltzman's constant, and N the number of lattice sites. 
Next, by putting our ansatz A(k, ui) and B(k, uj) into the set of moment equations, we get to lowest order in 



= i / n k - ffll (k) 
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a3(k) = 1 — ai(k) — a 2 (k) 



where the quasi-particle energy is given by 



with 



ff k ~ffo(k)=17i(k) 

n k « n (k) = n 2 (k) 
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where the f2j(k)'s are the ones which enter in the moment equations in the normal state (Eqs. ( [l7| ) and Ref. M ). 
However, due to the fact that the narrowing factor, B, is a sum over three frequencies of our problem (see Eq. (|14|)) 
then all three parameters (p,B and A(T)) are mixed together. 
From Eqs. 

we see that the quasi-particle energy, E k , has a gap which we find to be: 



A(T,k) = A (T) 



» 2 (k) -fli(k) 
f2 2 (k)-ai(k) 



(18) 



From Eq. (jl|) we conclude that our gap equation, which is a manifestation of singularities in the density of states, 
is k-dependent. Going back to our local Hamiltonian (see Eq. (0)), which in reciprocal space is a constant, then at 
the mean field level we should obtain a pure s-wawe, while going beyond mean field approximation we have been able 
to derive a k-dependent gap. In our case, we have been able to modify the BCS results, from pure s-wave to a wave 
vector dependent energy gap. 

Our Ansatz is based on the assumption that the role of correlations is mainly taken into account in the diagonal 
one-particle spectral function. This implies that the off-diagonal order parameter has been taken to as the BCS one, 
i.e., to m^^(k). In a previous work, jl0| it is found that both the diagonal and off-diagonal spectral functions have 
four peaks, symmetric in pairs, for U/t = —4. However, in the diagonal spectral function the fourth peak has a small 
weight and it is neglected here. If we want to consider two other peaks for the off-diagonal spectral function we have 
to include more moments, which is beyond the present work. Also, we have neglected life-time effects in Eq. (|l5|), 
which would require the evaluation of higher moments. We leave this out from the present calculation. 

In Fig. 1 we show the dependence of the reduced gap as function of reduced temperature. For comparison, we have 
also included the BCS gap. The critical temperature, T c /t, has the value 0.80 and the gap at zero temperature is 
0.66£. Then, for the ratio A(0)/T c we find 0.825. We should compare this with the BCS universal value, 1.76, for the 
same parameters [pT|- p"2] ]. In Fig. 2 we show the energy spectra for the three bands involved in the diagonal spectral 
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density. We observe that the third band, f2k, is almost always on top of the superconducting bands, i.e., ±i?k- In the 
inset of Fig. 2 we include a blowup of close to its minimun, in order to calculate the gap. In Fig. 3 we present 
the k dependent factors, 71 (k) and 72 (k), 



7i (k) 



(» 2 (k)-^(k))Ag(r) 
(n 2 (k)- ai (k))a 2 (k) 



72 (k) 



(n 2 (k)-n 1 (k))n 1 (k)Ag 

(fi a (k)-oi(k))o3(k) 



(T) 



(19) 



for T/t = 0.001, U/t = —4.0 and p = 0.1. These factors, which represent a correction to our frequencies, -ffk and 
f2k, are zero at T c . This implies that they do not change the value of T c . However, for T « they reach up to a 
.25. This contribution may look small but we have to remember that our equations are highly non-linear and a small 
change could produce a large change in the variables of interest, i.e., A(0). In consequence, our results support the 
calculations of Refs. |ll|- ]l2] in the sense that the critical temperature is renormalized by more than a factor of two 
with respect to the BCS value. On the other hand, we argue that the value of the gap at zero temperature is going to 
be modified by the inclusion of the 7(k)'s, which are going to increase the value of A(0). Due to these arguments, our 
calculation is of a perturbative character in A^ (T) . In the calculation of Pedersen et al |l(| the authors find a critical 
temperature of T c = 0.19 by using a different moment approach where the spectral functions have four peaks. Our 
calculations for the spectral weight supports the view of ref |n| because there is an additional band which splits off 
from 0:3 (k). []l3| The appearence of a fourth band will modify the value of T c and also the value of the gap. However, 
the numerical Tc/t in the present work is different than the one given in ref. JTo[ ] . 

By using the moment approach in the presence of off-diagonal long range order, i.e., A D (T) ^ 0, to the negative 
Hubbard model, we have worked out the attractive Hubbard model in 2-D with an ansatz where the diagonal spectral 
function has three peaks and the off-diagonal one has two. The physical meaning of each of these peaks has been 
discussed. We have solved our equations in lowest order of A Q (T), an approximation which obviously is valid close to 
T c but it fails for zero temperature. We have seen that the effect of the third band is to renormalize the order parameter 
producing an energy gap which is k-dependent. In consequence, in our approach the order parameter, A (T), and 
the energy gap, A(T, k), are not the same. A similar discussion was presented by Randeria, Duan and Shieh [ jl4| . 
The implications of this point is that we can shed light on the experiments carried out in the high-temperature 
superconductor materials since the experimentalists are talking about the symmetry of the order parameter and from 
the theoretical side we should refer to the energy gap. These two concepts are valid in mean field theory, i.e., in the 
BCS approximation. When we go beyond BCS as in the present work, it is difficult to asses the symmetry of the order 
parameter from tunelling experiments only. Then, from the experimental side it is difficult to find the symmetry of 
the order parameter unless the high temperature superconductors obey mean field equations. In short, the messages 
of this paper are 

• 1- Due to correlations, the order parameter and the gap are two different quantities. In our case, from a local 
pairing attraction, i.e., pure s-wave in mean field, we get a k gap function. 

• 2- Our calculations are of a perturbative character and they agree with Rcfs. |llj- Work is in progress fi"3j ] 
to include additional bands in the diagonal spectral function. This will be published elsewhere. 
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Figures. 

Fig. 1.- The energy reduced gap, A(T)/A(0), as function of reduced temperature, T/T c , for U/t = —4.0, p = 0.1. See 
Fig. 2 (inset) for the definition of A(T). 

Fig. 2.- The energy bands, ±£"k and Ok along the diagonal of the Brillouin zone. For comparison, we have included 
the BCS energy bands. In the inset we have made a blow up of close to the minimun of the band. This allows 
us to find A(T). Same parameters as in Fig. 1 with T/t = 0.001. 

Fig. 3.- We plot 71,72 vs k, along the diagonal of the Brillouin zone for the same parameters in Fig. 2. The parameters 
7's are defined in Eq. (19). 
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